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Yuki IzumidgEI and Koji Okuda 

Department of Information Sciences, Ochanomizu University, Tokyo 112-8610, Japan 
Division of Physics, Hokkaido University, Sapporo 060-0810, Japan 

We formulate an endoreversible finite-time Carnot cycle model based on the assumptions of local 
equilibrium and constant energy flux, where the efficiency and the power are expressed in terms of 
the thermodynamic variables of the working substance. By analyzing the entropy production rate 
caused by the heat transfer in each isothermal process during the cycle, and using an endoreversible 
condition applied to the linear response regime, we identify the thermodynamic flux and force of 
the present system and obtain a linear relation that connects them. We calculate the efficiency at 
maximum power in the linear response regime by using the linear relation, which agrees with the 
Curzon-Ahlborn efficiency known as the upper bound in this regime. This reason is also elucidated 
by rewriting our model into the form of the Onsager relations, where our model turns out to satisfy 
the tight-coupling condition leading to the Curzon-Ahlborn efficiency. 

PACS numbers: 05.70.Ln 


I. INTRODUCTION 


The physics of heat engines originates from the Carnot’s great discovery of the fundamental upper bound of the 
thermodynamic efficiency 77 of heat engines working between two heat reservoirs with temperatures and Tf- 

{Tl>Tf) dHi: 

< 1 - 7 ^ = ? 7 C (Carnot efficiency), (1) 

where the equality holds for an infinitely slow process (quasistatic limit) with zero dissipation realized in, e.g., the 
Carnot cycle. Indeed, this discovery may be regarded as the origin of thermodynamics itself. However, the quasistatic 
limit is an ideal case, and the thermodynamic processes observed in daily life occur at finite rates. Remembering 
that we always demand power for our use of electric devices, which may originally be generated from power plants 
converting heat flux into electric power, we require the physics of powerful heat engines free from the limitation of 
the equilibrium thermodynamics. This deep understanding of powerful heat engines is becoming more important due 
to the worldwide energy crisis and climate change. 

The physics of heat engines maximizing the power rather than the efficiency was developed in a classical paper by 
Curzon and Ahlborn [^. They showed that, under an assumption of the endoreversible condition and the Fourier law 
of heat conduction between the working substance and the heat reservoir, the efficiency at maximum power rj* of a 
finite-time Carnot cycle is given by the following Curzon-Ahlborn (CA) efficiency: 


fr^ AT^ 

= = '2f^ ~ + ■ ■ ■ ’ 

where we have defined the temperature difference between the heat reservoirs as AT^ = and the averaged 

temperatures between the heat reservoirs as ^ ^ ^ . We note that the same formula was derived more previously 

in, in 1950s, and it is claimed that its origin even goes back to a paper in 1920s according to a recent 

work 0. Because this CA efficiency displays a similar simplicity to the Carnot efhciency, it led to the development 
of a new discipline of finite-time thermodynamics that aims to account for the efficiency of actual power plants and 
thermal devices [iSl. The key to the derivation of the CA efficiency is the phenomenological assumption of the 
endoreversible condition (as named by Rubin Q later), which means that the irreversibility occurs only by the heat 
transfer process between the working substance and the heat reservoir, and that the state of the working substance 
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is internally reversible whose entropy change along the cycle is expressed by a Clausius-like equality (see Eq. (fT3l) b 
Under this condition, the efficiency of the finite-time Carnot cycle is given by the Carnot efficiency using the ratio of 
the temperatures of the working substance (see Eq. (114^ 1. and the efficiency at maximum power rj* is expressed by 
using the square root of the temperature ratio T^/T^ as in Eq. @ i- 

Despite its importance, even until recently, there has been no argument showing whether the CA efficiency is 
universal as rj* from the viewpoint of fundamental physics. The role of the CA efficiency has become increasingly 
important after Van den Broeck [T^ proved that the CA efficiency is the upper bound of rC in the linear response 
regime by using the Onsager relations of the linear irreversible thermodynamics framework [l6|-[l8j: 


r,‘<^=VCA + 0 (at-") , 


( 3 ) 


where the bound is realized under the tight-coupling (no heat-leakage) condition [H- Because the Carnot efficiency 
Eq. ([T]) can be approximated as rye — in the linear response regime, it is also equivalent to say that rj* in the 
linear response regime is bounded from above by one half of the Carnot efficiency as 77 * < ^. As reviewed in di¬ 
ll^, since the paper by Van den Broeck [l^ . various studies on finite-time heat engines have been conducted, which 
include linear response [23 - [30l| , nonlinearres 2 onse , stochas tic |4ll-|43 , quantum , thermoelectric 


photoelectric [^, molecular dynamics |52l - l55l| . and experimental |56l-l58l| studies. 

In , Van den Broeck established a view that the process of heat energy conversion into work in the linear response 
regime is ruled by a cross effect based on the Onsager relations: 


Ji = LiiXi + L12X2, 
J2 = L21X1 + L22X2, 


( 4 ) 

( 5 ) 
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an identical formulation has also 


where Xi is an “external” thermodynamic force and X 2 is a “thermal” thermodynamic force that is proportional to 
AT^, Ji and J 2 are their conjugate thermodynamic fluxes, and L^-’s are the Onsager coefficients with reciprocity 
Li 2 = L 21 (see Sec. II C for details). While this viewpoint is familiar in steady-state heat energy conversion, such as 
in thermoelectric devices usually analyzed with the Onsager relations [J, 1^ 

been established even for cyclic heat engines such as a finite-time Carnot cycle 

Despite these successes, these theories of finite-time heat engines may still be abstract compared to the theory of 
the quasistatic heat engine. One reason could be that a general state of a finite-time heat engine cannot be drawn on 
a thermodynamic plane, thus a clear picture is lacking, unlike the quasistatic cycle with well-defined thermodynamic 
variables of the working substance such as the pressure, volume, and so on. Even in a finite-time cycle, however, it 
may still be possible to assume that the state of the working substance is specified by a unique combination of the 
thermodynamic variables at any instant along the cycle, and hence the working substance and the heat reservoirs 
are in a local equilibrium, but not in a global equilibrium similar to the quasistatic cycle. According to this local 
equilibrium assumption, we can draw the finite-time cycle on the thermodynamic plane as well as the quasistatic 
cycle, and can also assume that the fundamental thermodynamic relations hold between the thermodynamic variables 
of the working substance during a finite-time cycle. In fact, Rubin introduced this local equilibrium thermodynamic 
description to the endoreversible cycle very previously Q (see also [1^), where the endoreversible condition is expressed 
by the entropy change of the working substance during one cycle. However, this local equilibrium assumption has not 
fully been taken into account in the recent literatures, and we hence are naturally motivated to elucidate how these 
endoreversible heat engine models based on the local equilibrium assumption are unified with the more recent linear 
irreversible thermodynamic description using the Onsager relations Eqs. ([4]) and ((5|) . 

In the present study, we formulate an endoreversible finite-time Carnot cycle model based on the assumptions of 
local equilibrium and constant energy flux. In our framework, the power and the efficiency can be expressed in terms 
of the thermodynamic variables of the working substance. Erom the analysis of the entropy production rate caused 
by the heat transfer in each isothermal process during the cycle, we identify the thermodynamic flux and force in each 
isothermal process, where the flux and force are assumed to be related by the Fourier law. We also find that, due 
to the endoreversible condition applied to the linear response regime, these thermodynamic forces in the isothermal 
processes are not independent, and we can identify the reduced thermodynamic force, whose conjugate thermodynamic 
flux turns out be an averaged heat flux. We then obtain a linear relation that connects these thermodynamic flux 
and force. From the calculation of the efficiency at maximum power by using this linear relation, we obtain the CA 
efficiency as 77* of the present model. Then, suitably changing the variables, we elucidate that the linear relation 
in our framework can be rewritten into the form of the Onsager relations, from which we can directly confirm the 
tight-coupling condition leading to the CA efficiency in the framework of [1^ . Therefore our work establishes a precise 
connection between the finite-time thermodynamic approach to heat engines by Curzon and Ahlborn and the linear 
irreversible thermodynamic approach based on the Onsager relations [l5l | via the local equilibrium assumption that 
gives us a more intuitive picture of the finite-time heat engines. 
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The present paper is organized as follows. In Sec. II A, we introduce our model based on the assumptions of local 
equilibrium and constant energy flux. In Sec. II B, by analyzing the entropy production rate of our heat engine with 
the aid of the endoreversible condition applied to the linear response regime, we naturally identify the thermodynamic 
flux and force and a linear relation that connects them. We then calculate the efficiency at maximum power in the 
linear response regime by using this linear relation. In Sec. II C, we elucidate the relationship between the linear 
relation obtained in Sec. II B and the Onsager relations Eqs. m and explicitly showing that our model surely 
satisfies the tight-coupling (no heat-leakage) condition leading to the CA efficiency. In Sec. Ill, we discuss a few 
aspects related to our formulation in Sec. II, and summarize our study. 


II. MODEL AND RESULTS 


A. 


Local equilibrium thermodynamic formulation of the endoreversible finite-time Carnot cycle 


Our heat engine model consists of the working substance, the hot heat reservoir with temperature and the cold 
heat reservoir with temperature T^. We assume that the working substance is always in a local equilibrium state 
specified by a unique combination of the well-defined thermodynamic variables, and the heat reservoirs are also in a 
local equilibrium state. Denoting the internal energy and the entropy of the heat reservoir by C//*' and {i = h,c), 

1 dS^ 

respectively, the temperature is defined by ^ . We also denote the internal energy and entropy of the working 

substance by U and S, respectively. Hereafter, we use the suffix i to denote the thermodynamic variable of the 
working substance when it contacts with the heat reservoir with the temperature We then obtain the first law of 
thermodynamics (energy-conservation law) by using these thermodynamic variables as @ 


dC/f _ dU, ^ dVi 
dt dt ^ dt ' 


( 6 ) 


where we define ^ 


Q Q -p /DO 

= and through the fundamental thermodynamic relation 

dU,{S,, V,) = ^dS, + = T,dS, - P^dV,, 

aSi dVi 


( 7 ) 


with Ti, Pi, and Vi being the temperature, pressure, and volume of the working substance, respectively. Eq. states 

dU^ 

that the heat flux from the heat reservoir —which is the internal-energy change rate of the heat reservoir, is 

decomposed into the internal-energy change rate of the working substance and the instantaneous power output 
p m 

Our heat engine experiences a thermodynamic cycle that consists of (I) an isothermal expansion process in contact 
with the hot heat reservoir with the temperature T^, (II) an adiabatic expansion process, (III) an isothermal com¬ 
pression process in contact with the cold heat reservoir with the temperature and (IV) an adiabatic compression 
process (see Eig. [T](a)). We assume that we can take the durations of the adiabatic processes to be sufficiently short 
compared to the ones of the isothermal processes in a finite-time cycle, and that the thermodynamic states of the 
working substance move along the quasistatic adiabatic curves in the thermodynamic plane, while the thermodynamic 
states of the working substance during the isothermal processes do not agree with the quasistatic isothermal curves. 
This assumption is quite natural because relaxation to the equilibrium state during the adiabatic process is occurred 
by internal processes inside the working substance itself, which is much faster than the speed of relaxation to the 
(global) equilibrium during the isothermal process determined by interactions with the heat reservoir. Erom this, our 
heat engine can be regarded to run one cycle in the cycle time tcyc = th P tc, where we denote by ti the duration 
of the isothermal process in contact with the heat reservoir at Additionally, in our formulation, we also assume 
that the energy flux corresponding to each term in Eq. is constant at any instant along each isothermal process, 
and the temperature Ti does not change during the isothermal process, where these assumptions are also adopted in 
the original CA model (see Sec. Ill for a discussion on this assumption). Using Eq. (O, we can rewrite Eq. ® as 


dUf ^ dS, 
dt * dt 

The heat from the heat reservoir during the isothermal process Qi is calculated by using Eq. m as 


( 8 ) 


Qh = 

Qc = 


I 

L 


th jtjR 

dt = ThASh, 


dt 

dU^ 


dt = TrASr 


(9) 


dt 


( 10 ) 
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FIG. 1: (a) Temperature-entropy (T-S) diagram of an endoreversible finite-time Carnot cycle based on the assumptions of 
local equilibrium and constant energy flux, where T and S are well-defined (the inner cycle depicted with the thin solid curve) 
on a thermodynamic plane even far from the quasistatic limit (the large outer cycle depicted with the bold solid curve). The 
cycle consists of (I) isothermal expansion process in contact with the hot heat reservoir with the temperature T^, (II) adiabatic 
expansion process, (III) isothermal compression process in contact with the cold heat reservoir with the temperature T^, and 
(IV) adiabatic compression process. In the inner finite-time cycle, the temperature of the working substance Ti during the 
isothermal process is assumed to be constant at any instant, but does not agree with tP'. The area of the cycle represents the 
work output W during one cycle; W = ATLAS' for the quaistatic cycle and W = AT AS for the finite-time cycle, (b) Schematic 
illustration of the heat energy conversion using the endoreversible heat engine model, where the working substance specified by 
a unique combination of the thermodynamic variables and the heat reservoirs are assumed to be in a local equilibrium state. 
The thermodynamic fluxes and forces in the figure are defined in Eqs. (HI), dMl) and dig). 


where we defined the entropy change of the working substance during the isothermal process as 

ASh = / —^dt, ASc = / —^dt. II) 

Jo dt ^ it, dt ^ ’ 

Because we require that the cycle is closed after one cycle and we also assume that the adiabatic processes are regarded 
as quasistatic processes, the following relations should hold: 


ASh = -ASc = AS. 

We note that the endoreversibility condition [ 1 , 


Qh Qc 


( 12 ) 


(13) 


automatically holds from Eqs. (0), dini) and m- This condition manifests that the entropy production is occurred 
only by the heat transfer between the heat reservoirs and the working substance, and the internal state of the working 
substance is, as assumed above, in a local equilibrium state, which implies that the entropy change of the working 
substance along the cycle is expressed by the Clausius-like equality using the temperature of the working substance. 
Eq. (fT51) also implies that the efficiency of this type of the heat engines is given by the “endoreversible Carnot 
efficiency” using the temperatures of the working substance as 


V = 


Qh Th 


(14) 


Although this differs from the usual Carnot efficiency Eq. ([T|) that uses the temperatures of the heat reservoirs, it 
is remarkable that the efficiency is still expressed by the thermodynamic variables of the working substance. In this 
way, the endoreversibility condition Eq. (1131) largely constrains the behavior of the heat engines. The power output 
W of the heat engine is also expressed by using the entropy change as 

Qh + Qc at AS 

f f f ^ 

^cyc ^cyc ^cyc 


( 15 ) 
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where AT = Th — Tc is the temperature difference of the working substance between the isothermal processes and 
we used the first law of thermodynamics W = Qh + Qc for one cycle. Hereafter we denote by the dot the quantity 
divided by the cycle period. Then the rightmost expression in Eq. (TTKll using the “endoreversible work” W = AT AS 
written by the thermodynamic variables of the working substance like the quasistatic work W = ATLAS' may be 
regarded as “endoreversible power.” The power Eq. (flSl) and the efficiency Eq. (HI, which are defined by using the 
thermodynamic variables of the working substance in this way, are a characteristic of our local equilibrium description 
of the endoreversible heat engine model. Especially, the “endoreversible power” expression has never been used in 
a context of determining the efficiency at maximum power so far. We note that these expressions also give us an 
intuitive picture underlying heat-energy conversions in a finite-time Carnot cycle, because the efficiency and the work 
output are still given by the quasistatic-like expressions. 


B. Efficiency at maximum power in the linear response regime 


In this subsection, we consider the efficiency at maximum power of our heat engine in the linear response regime 
AJ'R based on the local equilibrium assumption introduced in Sec. II A. 

Our analysis using the linear irreversible thermodynamics begins from the entropy production rate. Because we 
assume that the heat reservoirs and the working substance are always in a local equilibrium state with the well-defined 
entropies, the entropy production rate ^ of the total system (the working substance and the heat reservoirs) at any 
instant along the isothermal process is expressed by the sum of the entropy change rates of these partial systems: 


dt dt dt 


dS^ dV, 
dUi dt dVi dt 


I 1 


dt dt dUf- dt 

where we used Eq. ®. Then, the entropy production rate for one cycle a is written as 


dt ’ 


1 


tc 


ds = 


I 


t, 


eye Jo 


dt t, 


eye JO 


= E 


I 


1 


1 

n 

AUf^ 


dU^ , 1 

—7^dt-\ -/ 

dt tcyc Jth 


eye / ^ ^ 

W~Tc 


Tj tcyc 

where we defined the internal energy change of the heat reservoir during the isothermal process as 


(16) 


dUf 

dt 


dt 


(17) 


AU^ = 


n dt ’ 


AU^ = 


Jth 


dt 


-dt. 


(18) 


From Eq. (HZl, we naturally define the thermodynamic force as the (inverse) temperature difference between the 
working substance and the heat reservoir, and we define the conjugate thermodynamic flux as the heat flux from the 
heat reservoir (see Fig.[T](b)): 


Jq. 


AUf 


Qi Y 

T-, -^Ti 

^cve 


1 1 


(19) 


Because the energy flux and the temperature of the working substance is assumed to be constant along each isothermal 
process, we obtain 


AU^ dU^ ^ dS, 
ti dt ^ dt ^ 


( 20 ) 


by using Eq. ®. Then the thermodynamic flux Jq. can also be expressed by using the time derivative of the 
thermodynamic variable of the working substance as 


jQd = ( 21 ) 

where we denote by the ratio of the duration of each isothermal process ti to the cycle time tcyc- To proceed 
further, we need a relation that connects Jq^ and Xt^, in addition to the local equilibrium thermodynamic formulation 
mentioned in Sec. II A. Because we are adopting the local equilibrium thermodynamic assumption, it is also quite 
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natural to assume that the heat flows in proportion to the temperature difference (the Fourier law) in the same way 
as the original CA model 


Q, = k,{T^ - T,)U, (22) 

where we denote by Ki the thermal conductance between the heat reservoir with the temperature Th and the working 
substance. Using Eq. (|2^ , we then obtain the following relationship between Jq^ and Xxi ■ 


Jq. 


^ = aMTl" - Ti) = a,K,T^T,XT^. 

^cyc 


(23) 


In the following, we consider the linear response regime AT^ —>■ 0. By substituting Eq. (|2^ with Citcyc instead 
of ti into the endoreversibility condition Eq. m, and expanding that with respect to AT and AT^, we can write 
Eq. (fOl) as 


/T^ \ ( AT AT^\ 

{ahKh + acKc) (—-!)+ {ahKh - a^Kc) ^ 1 


+ 0(AT^AT,AT2) =0, 


(24) 


where T = is the averaged temperature of the working substance between the isothermal processes. Then, by 

substituting T = +c^'^AT^ + c^^^AT into Eq. we obtain the expansion coefficients order by order as follows: 


c 


( 0 ) _ 


c 


( 1 ) 

R 


0>c^c 

^{ahKh + acKc) ’ 


c 


(1) = 


2{ahKh + OcKc) ’ 


Then, Th = T + and Tc = T — are expressed as 


T = — Qc^c 

2,(^(lh^h H“ ^c^c) 



Ojc^c 

‘ZijXh^h 0>cf^c) 
Clc^c 

‘^{p-h^h H“ ^c^c) 



(25) 

(26) 

(27) 

(28) 
(29) 


By using Eqs. ([25|)-(|27|1. we also find that Xt^ ’s are expressed in terms of AT and AT^ as follows: 


^Tc 


- Th 
T f - T, 


acKc AT^ — AT 

ahKh + acKc T^^ 

ahKh AT^ - AT 

ahKh + acKc T^^ 


acKc 


ahKh + acKc 
ahKh 


Xt, 


ahKh + acKc 


Xt, 


where we defined the “reduced thermodynamic force” Xt as 


Xt 


AT^ - AT 


(30) 

(31) 


(32) 


From Eqs. dsni) and dSU), we find that Xt^ and Xt^ are proportional with each other in the linear response regime. 
By using Eqs. o and (1311) and approximating Jg. in Eq. (1231) as Jq^ ~ UiKiT^ ^Ti, we can simplify the entropy 
production rate Eq. (ED by using Xt up to the quadratic order of Xt as 


I-eye 


ds = JQi^Ti — ^^aiKiT^ Xt^ 


where Jq is the averaged heat flux [39|, 1^, l60j defined as 


ahKh + acKc 


T^^Af = JqXt, 


T - JQh- JQc 

Jq = -^- 


- 1 Ax’. 

“t“ (Xq^c 


(33) 


(34) 
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Therefore the description of the present heat engine model is reduced to this linear relation Eq. (l34l) . In the lin¬ 
ear response regime, the endoreversible Carnot efficiency Eq. (imi and the endoreversible power Eq. (Ha are also 
approximated by using Xt as 


V 

W 


AT 

'jW 


AT^ 

J'R 


- T^Xt, 


- T ^ T - 

— ^Qh rpR I 

T Qj}iK}i -|- CLc^c 


T\ 


^Xt[ 


AT^ 



(35) 

(36) 


respectively, where VE is a quadratic function of the reduced thermodynamic force Xt. Because Xt is proportional to 
the temperature difference between the working substance and the heat reservoir as in Eqs. dsnu and m, it is natural 
that we can control the power by changing Xt- the quasistatic limit is realized under no temperature difference 
between them as Xt^ = 0 (AT'i® = AT^). Then the heat flux Jq^ also vanishes from the linear Fourier law Eq. (l2^ 
via Eqs. dSOl) and (ED, and we obtain no power VE = 0 from Eq. (15^ as in the usual Carnot cycle. As Xt increases 
from X!^^, the heat flux Jq. becomes finite. Since the power VE is a quadratic function of Xt, it can take a maximum 
value at a certain point, which is determined by = 0 : 


A 


* _ 

T — 


A^R 




(37) 


In the quasistatic limit X!^^ = 0, we attain the Carnot efficiency ~ 77 c as the maximum efficiency from Eq. (1351) . 
At the maximum power, from Eqs. (1551) and ED, the efficiency 77 * is given by 


a^r 

2 TR ■ 


(38) 


This is the CA efficiency up to to first order of AT^ that corresponds to the equality in Eq. ([3|) . The maximum power 
IE* is also given as 


ahKhttcKc AT^^ 
ahKh + acKc 


(39) 


which depends on the thermal conductivity [^. Here we note that we obtain the CA efficiency in Eq. (l38l) without 
using the Onsager relations Eqs. (j4|) and ([D, but using the linear relation between the thermodynamic flux and force 
Eq. (IM)l with the aid of the endoreversible expression of the efficiency and power Eqs. (IM|) and (|M)) . In the next 
Sec. II C, we consider this connection. 


C. Formulation of the endoreversible flnite-time Carnot cycle model using Onsager relations 

As we have shown in Sec. II B, the efficiency at maximum power 77 * in Eq. (1551) . which is based on the local 
equilibrium thermodynamic formulation using the linear relation Eq. (|34l) . is the upper bound in Eq. ©, while the 
inequality in Eq. ([D comes from the formulation based on the Onsager relations [1^. Therefore, we elucidate the 
relationship between these formulations in this subsection. 

First, we briefly review the derivation of the inequality for the efficiency at maximum power in Eq. mM- Denoting 
an external force and its conjugate variable by F and x, respectively, we can generally express the power of the heat 
engine IE as IE = —Fx. Then the entropy production rate of the total system & is decomposed into the sum of the 
entropy increase rate of each heat reservoir because the state of the working substance should return to the original 
state after one cycle: 


Qh Qc 

Jf-Jf 



J.R ) yR 


Fx 

J'R 


AT^ ■ 

Qh = JlXi + J2A2. 


(40) 


Here, the thermodynamic fluxes Ji = x and J 2 = Qh, and their conjugate thermodynamic forces Xi = ^ and 
X 2 = related through the Onsager relations Eqs. dD and ([5]). We note that, in this case, we do not necessarily 

assume that the working substance along the cycle is expressed in terms of the well-defined thermodynamic variables, 











in contrast to our formulation in Sec II B. Using these thermodynamic fluxes and forces, the power and the efficiency 
are given as 


w = -JiXiT^, 

(41) 

W JiAiTR 

(42) 

'^~Qh~ J 2 • 


With these expressions as well as the Onsager relations Eqs. ([4]) and (O, we find that the maximum power is realized 
at Xi = — from = 0. Its efficiency rj* is given as 


^ 2 - g2 2T^ ’ 

which is a monotonically increasing function of jgl, where the coupling strength q is defined by 


(43) 


9 = 


Li2 

'JLi\L22 


(44) 


From the non-negativity of the entropy production rate a = JiXi + J 2 X 2 , the Onsager coefficients L^’s should satisfy 
^11 > 0, L 22 > 0, and L 11 L 22 — L 12 > 0, and they impose the following constraint on q\ 


\q\ < 1, (45) 

where the equality is known as the tight-coupling (no heat-leakage) condition [l^, [l^. Under this tight-coupling 
condition, ry* in Eq. (14311 attains the upper bound given by the CA efficiency rjcA as in Eq. ©• An essential point 
of the derivation of the formula Eq. (j43l) is that the non-zero cross-coefficient L 12 plays an important role in iq* in 
Eq. (031), which is clear from the definition Eq. 04|). 

Returning to our original problem, from Eq. (1^51) . we formally obtain the following “Onsager coefficients” under 
our choice of the thermodynamic fluxes Jq^ and forces Xt^ '■ 


f 0 \ 

0 


(46) 


where there are no nondiagonal elements. This contrasts to the formulation using Eqs. o and where the cross¬ 
terms play an important role in the heat-flux conversion into power [H. Eq. (|46ll is natural if the entropy production 
originating from the heat transfer between the working substance and the heat reservoir in each isothermal process 
is independent of each other. However, as seen from Eqs. (0Ql) and m, Xxi’s are not independent of each other, 
unlike Aii’s in Eq. (IMl) . but At^’s and Xi’s should be related with each other by a variable change. Moreover, how 
the external force F in the definition of Xi = F/T^ against the heat engine is related to the thermodynamic 
variable of the working substance is not obvious in our local equilibrium description of the finite-time Carnot cycle 
where the efficiency and the power is expressed by them as in Eqs. da and a. To elucidate these points, we restate 
our expression of the entropy production rate Eq. a with At’/s using independent thermodynamic forces as (see 
Fig. [11(b)) 


1 1 \ AC/^ 


J'. 


tc. 


/ j_^ 

W T^J tcyc 


AT AS 

^jQ^Y^n + JsYT. 

C ^cyc 


(47) 


We can make this restatement by using the endoreversibility condition Eq. a, the first law of the thermodynamics 
bE = Qh + Qc for one cycle, and Eq. a, where we defined the heat flux from the hot heat reservoir as a new 
thermodynamic flux and its conjugate new thermodynamic force as 


JQh 


Qh 

f eye 


auTh 


dt 


Y-rn = 


TR 


I 

w 


(48) 


In addition, we defined the entropy flux as another new thermodynamic flux and its conjugate new thermodynamic 
force as 
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In this way, all the thermodynamic fluxes and forces are expressed in terms of the combination of the thermodynamic 
variables of the working substance and the heat reservoirs owing to the local equilibrium assumption. From Eq. (ggi, 
in particular, the new thermodynamic force Yt is proportional to the temperature difference of the working substance 
between the isothermal processes (Fig. [T] (b)), which has never been proposed so far. The external force F in the 
linear irreversible thermodynamics framework in Eq. m by which the power is maximized is abstract in the case of 
the finite-time Carnot cycle. But our definition Yt in Eq. (|4^ gives us a more intuitive picture of “force” applied to 
the finite-time Carnot cycle, which can be used as a parameter for maximization of the power. This is also consistent 
with the Curzon and Ahlborn’s original idea that reminded them of the maximization of the power in their finite-time 
Carnot cycle [^: the power becomes zero in the quasistatic limit AT = AT^ and also in the other extreme situation 
AT = 0 where the heat flux from the hot heat reservoir is just transferred to the cold heat reservoir without generating 
any work. Then the power must take a maximum with an intermediate AT between these two zeros. 

In the linear response regime, these new thermodynamic fluxes and forces are approximated as 






dt 


Yj'R 


AT^ 


Js 


AS dSh 

tcyc dt 


Yt 


AT 

yR ■ 


(50) 


We note that Yt (Ttr) corresponds to Xi (A 2 ), and Js {Jqh) does to Ji (J 2 ) in the definition of the general 
formulation in Eq. (jM]), respectively. From Eq. ([50]), we immediately notice that Jq^ and Js are in proportion to 
each other: 


Jq. = T^Js- (51) 

In fact, the proportionality between the two thermodynamic fluxes in Eq. (1511) indirectly implies the tight-coupling 
condition of this system |g| = 1 [s^, because we can easily show from Eqs. (|4|) and ([5]) the relation J 2 = + 

T22(1 — q^)X 2 between the two thermodynamic fluxes. However, to understand this relationship more directly and 
precisely, we express the present system by the following Onsager relations using the new thermodynamic fluxes and 
forces: 


Js = LttYt ~\~ Ltt'^Yt^i (52) 

= Lt^t^t T (53) 

To obtain the new Onsager coefficients from the previous coefficients in Eq. (l46|) . we relate the thermodynamic forces 
XTi {i = h,c) and Y^ {m = T, T^) by using Eqs. (1501) . (I5T|) . and (150)) as follows: 

[/ QcKc \ ^ _ I _ f acKc acKcT^ ^ ^ (bdl 

atiKh + acKc \—o.hK.h) {ahUh + aci^c)T^ \—ahKh —ahKhT^J \ YtrJ 
Rewriting Eq. (1541) as XTi = PirnYm hr Einstein notation, we obtain the new Onsager matrix 


_ ahKhacKc /I T^\ 
ahKh + acKc T^^J 


(55) 


from the relation Lmn = that conserves the entropy production rate as tr = LijXTiXT^ = LmnYmYn- 

Alternatively, we can directly obtain the Onsager coefficients Ltrt s-nd LtrtR from the expression of 

T rpR2 ^ tthKhClcKc rpR2 Y O-h^hdcKc rpR,^ . dh^httcl^c 

jQh='^hHhl ATfe = -;- 1 Xt = -;- 1 FT H-;-J Fj^R , (5Dj 

U/jR/i -(- Qjcl^c T 

which is obtained using Eqs. (1301) and (I54L and we can also obtain Ltt and LttR from Eq. m- 

Js = ^ = ahKhT^XT, = .^h^lh^sH^T^XT = ^hUh^^YT + .^hUh^^T^YTR. (57) 

T^ ahKh + acKc ahKh + OcKc ahKh + acKc 

From Eq. (1551) . it is straightforward to conhrm that the Onsager reciprocity and the tight-coupling (no heat-leakage) 
condition are fulfilled: 


LttR — Ltrti Q — 


LttR 


= I. 


LtRtR ^'T'T 


(58) 
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This implies that our model attains the upper bound (the CA efficiency) in Eq. ([3]) corresponding to |g| = 1 in 
Eq. (1431) . from the viewpoint of the linear irreversible thermodynamics framework. For completeness, we explicitly 
confirm that the maximization of the power by Yt leads to the CA efficiency by using the Onsager relations along with 
the general formulation introduced in the beginning of this subsection. Let us compare the two ways of maximization 
of the power by Xt as in Sec. II B and by Yr- As we have seen in Eq. (1541) . the reduced thermodynamic force Xt 
introduced in Eq. (|32]) is linearly connected by Yt and Ytr. By substituting Eq. (|54)l into Eq. we can express 
the power W in Eq. (1551) by using Yt and Ytr instead of using Xt as 


W = - 


Clc^c 


T^Y^ - 




(59) 


which can also be obtained from the expression W = —JsYtT^ in Eq. (HTl) and Eq. (1551) . We find that the quasistatic 
limit X^^ = 0 corresponds to = —T^Ytr from W = 0 by using Eq. (1551) . We also find that the maximum power 
point X^ = in Eq. (l37l) also corresponds to Y^ = —^Ytr from = 0 by using Eq. (l59l) . from which we can 
obtain the efficiency rj* = —Js^t'^^/'^Qh using Eqs. (1551) and (1551) as 

AT^ / 9\ 

^* = ^ = ^1ca + 0(at^ ), (60) 

Therefore, for our endoreversible heat engine model based on the local equilibrium assumption, we conclude that 
the efficiency at maximum power attains the upper bound in Eq. m from a viewpoint of the linear irreversible 
thermodynamics framework using the Onsager relations [H. We note that the above derivation is quite general 
because it does not rely on any particular working substance or thermal conductance. 


III. DISCUSSION AND SUMMARY 


We discuss a few aspects related to our formulation in Sec. II. 

First, we note that the endoreversible power as determined by the product of the heat flux and the endoreversible 
efficiency in Eq. (1551) in Sec. II B should also be determined by the first law of the thermodynamics W = Jq^ + Jq^ 
for one cycle. To this end, we need to consider Jq- with higher-order corrections of AT and AT^, while we have 
considered only the lowest order in Sec. II B. We utilize the expressions Jq^ = T^Js and Jq^ = —TcJs, which are 
obtained from Eqs. (ITUl) . (1T51) . (1T51) . and (1551) . Then, by using Eqs. (155|) . (1551) . and (1^ . we directly obtain the explicit 
form of Jq^ with the nonlinear terms as 


jQh 


JQc 


- 1 Xt -f 

O-hKh + OjcHc 

ahUh -b QcKc 


dhKhClcKc _ O-hl^hO'c^c 

2{ahKh + acKc) {ahKh + acKc)"^ 

|_ O'hl^hO'cl^c T^AT^X _ 

2{ahKh + acKc) (ahKh + Oc^c) 


Y'2 

■^T-> 

(61) 

T^^X^, 

(62) 


from which we can confirm W = Jq^ + Jq^ by using Eq. (1551) . 

Second, we point out that an extension of our formulation introduced in Sec. II A to the general case where the 
energy flux may not be constant along the cycle may be an interesting challenge. This would be more clarified by 
considering under what conditions the assumption of the constant energy flux strictly holds in a specific example: let 
us consider the following first law of thermodynamics Eq. m for the 3-dimensional ideal gas during an isothermal 
process: 



dT, _ 
^~dt ~ 




Td 


Nk^T, dV, 


(63) 


where we used the internal energy 14 = ^Nk-^Ti and the equation of state Pi = for the ideal gas, where 

we denote by N and /cr the particle number and the Boltzmann constant, respectively. For each term in Eq. (1631) 
to be constant, Ti = const solution to Eq. (1551) can be realized only under a specific protocol Vi{t) that satisfies 
= const. Therefore even in this simplest case of the ideal gas, our assumption of the constant energy flux 
introduced in Sec. II A (as also assumed in the original CA model) seems too restricted. In addition, although we 
assumed that the thermal conductance Kt is a constant in the present model (as also assumed in the original CA 
model), it can also depend on the thermodynamic variables such as Ti and 14 [27l . 1^ , and hence on time through 
them in general. Since actual heat engines may consist of complicated combinations of these factors, it is difficult to 
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specify generic conditions for arbitrary working substances and thermal conductances under which the assumption 
of the constant energy flux holds in a strict manner. Instead, for more natural description of the endoreversible 
finite-time Carnot cycle under more generic conditions, we are motivated to extend our formulation introduced in 
Sec. II A to non-constant energy flux cases that include the constant energy flux case as a special case. We are 
also in progress toward this goal from the linear irreversible thermodynamics point of view, which will be presented 
elsewhere [^. However, we stress that our present formulation of the constant energy flux case would serve as the 
“zeroth approximation” or a starting point for analyzing the performance of such actual heat engines. Indeed, this is 
also supported by the fact that observed data of the efficiency of actual heat engines are distributed around the CA 
efficiency [s^, which implies that they may be regarded as more or less variants of the constant energy flux case. 

In the present study, we formulated an endoreversible finite-time Carnot cycle model based on the assumptions of 
local equilibrium and constant energy flux. In our framework, the power and the efficiency are expressed in terms 
of the thermodynamic variables of the working substance. From the analysis of the entropy production rate caused 
by the heat transfer in each isothermal process, we identified the thermodynamic flux and force in each isothermal 
process, which are related by the Fourier law. By applying the endoreversible condition to the linear response regime, 
we found that those thermodynamic forces are not independent, and identified the reduced thermodynamic force and 
its conjugate thermodynamic flux, which are connected by the linear relation different from the Onsager relations. We 
calculated the efficiency at maximum power by using this linear relation, and obtained the Curzon-Ahlborn efficiency. 
We also elucidated that by suitable change of the variables, the linear relation between the thermodynamic flux and 
force in our framework can be rewritten into the form of the Onsager relations, where the novel thermodynamic 
force that is proportional to the temperature difference of the working substance between the isothermal processes 
is introduced. Then we directly confirmed that our model satisfies the tight-coupling condition that ensures the 
Curzon-Ahlborn efficiency as is the upper bound in the linear irreversible thermodynamics framework. We stress that 
our framework is quite universal because it only assumes that the working substance is in a local equilibrium state 
specified by a unique combination of thermodynamic variables at any instant along the cycle. We expect that our 
study unifies recent development of the theories of heat engines based on the universal nonequilibrium thermodynamics 
framework and the more phenomenological finite-time thermodynamics approach that was developed for application 
to real power plants and heat devices. 
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